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Abstract. Using high statistics samples of charged current interactions, MiniBooNE reports a model independent measure- 
ment of the single charged pion production to quasi-elastic cross section ratio on mineral oil without corrections for pion 
re-interactions in the target nucleus [1]. The result is provided as a function of neutrino energy in the range 0.4 GeV < E < 2.4 
GeV with 11% precision in the region of highest statistics. 



INTRODUCTION 

Future neutrino experiments will operate in the few-GeV neutrino energy range, where quasielastic scattering (CCQE) 
and pion production channels are dominant interactions. The charged current pion production interactions account for 
the biggest background for the neutrino oscillation disappearance measurements. Recently, three experiments (K2K 
[2], MiniBooNE [5, 6] and SciBooNE [3, 4]) have presented results for CC pion production cross section in the 
few-GeV neutrino energy region on nuclear targets. With high statistic samples it is possible to determine neutrino 
cross sections in this region with better precision than data from previous experiments [7, 8, 9] . Here, we present 
the ratio of cross sections of the charged current single pion production (CCl;r+, v^X jX^n^X') and quasielastic 
scattering {CCQE, v^n — > p) as obtained from the MiniBooNE experiment. In this measurement the largest sources 
of uncertainty coming from the neutrino flux determination largely cancels out. 

The Booster Neutrino Beam at Fermilab provides a neutrino source which is particularly well-suited to make this 
measurement; about 40% of neutrino interactions in MiniBooNE are expected to be CCQE and 24% CCn^ . The 
beam itself is composed of 93.6% with a mean energy of about 800 MeV and 5.9% (0.5%) Vfi (Vg) contamination 
[10]. 

The neutrinos are detected in the MiniBooNE detector [11], a 12.2 m diameter spherical tank filled with 818 tons 
of undoped mineral oil located 541 m downstream of the beryllium target. At the energies relevant to this analysis, 
the products of the interactions produce primarily Cerenkov light with a small fraction of scintillation light [11]. The 
light is detected by 1280 8-inch photomultiplier tubes (PMTs) which line the MiniBooNE inner tank. This inner tank 
region is optically isolated from a surrounding veto region, instrumented with 240 PMTs, that serves to reject incoming 
cosmic rays and partially contained neutrino interactions. 

Neutrino interactions within the detector are simulated with the v3 NUANCE event generator [12]. CCQE inter- 
actions on carbon are generated using the relativistic Fermi gas model [13] tuned to better describe the observed 
distribution of CCQE interactions in MiniBooNE [14]. Resonant CCln^ events are simulated using the Rein and Seh- 
gal (R-S) model [15], as implemented in NUANCE with an axial mass M]^ = 1.1 GeV . The angular distribution of the 
decaying pions in the center of mass of the recoiling resonance follows the helicity amplitudes of [16]. In MiniBooNE, 
87% of CCl7r+ production is predicted to occur via the A(1232) resonance, but 17 higher mass resonances and their 
interferences, as well as a nonresonant background that accounts for roughly 6% of CCln^ events, are also included 
in the model. Coherently produced CCl 7r+ events are generated using the R-S model [16] with the R-S absorptive fac- 
tor replaced by NUANCE'S pion absorption model and the overall cross section rescaled to reproduce MiniBooNE's 
recent measurement of neutral current coherent tt" production [17]. Coherent 7i:+ production is predicted to compose 
less than 6% of the MiniBooNE CCln^ sample due to the small coherent cross section [18, 19] and the dominance of 
the A++ resonance. A GEANT3-based detector model simulates the response of the detector to particles produced in 
these neutrino interactions. 
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ANALYSIS 



The CCQE and CCln^ reconstruction requires a detailed model of light production and propagation in the tank to 
predict the charge distribution for a given vertex and muon angle. The muon vertex, track angle, and energy, are found 
with a maximal likelihood fit, with the energy being determined from the total tank charge [20] . In the case of CCl7r+ 
events, only one track is reconstructed and assumed to be a jU^. The neutrino energy for both samples is reconstructed 
from the observed muon kinematics, treating the interaction as a 2-body collision and assuming that the target nucleon 
is at rest inside the nucleus: 

1 2mpE^,+m\-ml-ml 
£v=2 P— (1) 



nip — Efi cos 6 -^Eji — 



Here nip is the mass of the proton, is the mass of the muon, mi is the mass of the neutron in CCQE events and 
of the A(1232) in CCl;r+, 0^ is the reconstructed angle of the muon with respect to the beam axis (in the lab frame), 
and E^ is the reconstructed muon energy. 

The cross-section for a neutrino interaction is defined as: 

Nx 

= ^ J (2) 

^v^targets 

where N is the true number of events of type X(CC1 ;r+ or CCQE), Ntargets is the number of nuclear targets in the 
detector, and <I> is the flux of incident neutrinos. Because the incident neutrino flux and the number of nuclear targets 
are fixed properties of the beam and the detector (i.e. independent of the process we are considering), Ntargets and 
must be the same for both processes. Therefore we can write 



CTccpip{Ev) - (^ccqeiEv)- (3) 

With all the correction terms put together, the cross section ratio in each energy bin i is: 



(4) 



where the subscript ; runs over bins in true neutrino energy, subscript j indexes bins in reconstructed neutrino energy, 
Nx-cuts denotes the number of events passing cuts for X = CCln^, CCQE, f denotes the signal fraction, e denotes 
the cut efficiency, and f/ is a neutrino energy unsmearing matrix that acts on a reconstructed distribution to return the 
true distribution. 

To map reconstructed to true energy, we form a migration matrix Uij representing the number of MC events in 
bin / of reconstructed energy and bin j of true energy. We then normalize each reconstructed energy bin to unity 
to obtain an unsmearing matrix. This is equivalent to a Bayesian approach discussed in [21]; it differs from the 
standard matrix inversion method in that the resulting unsmearing matrix is biased by the MC distribution used to 
generate it. We account for this in our uncertainties by including a variation in the MC distribution used to generate 
the matrix. Because we have good data/MC agreement, this effect is small. The advantage of this method is that it 
avoids the problems of numerical instability and the magnification of statistical errors which occur in matrix inversion. 
This unsmearing procedure also proved insensitive to variations in neutrino energy reconstruction, confirming that it 
performs as intended. 



RESULTS 

For our primary measurement, we define CCl ;r+ signal as events with exactly one jj.^ and one 7r+ escaping the struck 
nucleus (which we call CCl7r+-like events) and CCQE signal as those with exactly one ji and no pions (CCQE- 
like events). Both event classes may include any number of protons or neutrons, but no other types of hadrons. The 
observed cross section ratio is then defined as the ratio of CCljC^-Mke to CCQE-]ike events and thus has not been 
corrected for re-interactions in the struck nucleus. 
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TABLE 1. The MiniBooNE measured CCl7i:+ to CCQE (as in Figure 1 but without the 
isoscalar correction) and CCl;r+-like to CCQE-like (Figure 1) cross section ratios on CH2 
including all sources of statistical and systematic uncertainty 



77 

(Itcv; 


cci/T fL,L,{j,t, yrav correciea^ 


cci/T — iiKe/L-L-{ii^ — iiKe \OoseT\Ga) 


0.45 ±0.05 


0.045 ±0.008 


0.036 ±0.005 


0.55 ±0.05 


0.130±0.018 


0.100±0.011 


0.65 ±0.05 


0.258 ±0.033 


0.191 ±0.019 


0.75 ±0.05 


0.381 ±0.047 


0.278 ±0.028 


0.85 ±0.05 


0.520 ±0.064 


0.371 ±0.040 


0.95 ±0.05 


0.656 ±0.082 


0.465 ±0.053 


1.05 ±0.05 


0.784 ±0.100 


0.551 ±0.066 


1.15±0.05 


0.855 ±0.1 14 


0.607 ±0.077 


1.25 ±0.05 


0.957 ±0.132 


0.677 ±0.091 


1.35±0.05 


0.985 ±0.141 


0.700 ±0.097 


1.5±0.1 


1.073±0.157 


0.777 ±0.109 


1.7±0.1 


1.233 ±0.207 


0.904 ±0.137 


2.1 ±0.3 


1.318±0.247 


1.022 ±0.161 



Figure 1 shows the observed CCl7r+-like to CCQE-like ratio extracted from the MiniBooNE data, including 
statistical and systematic uncertainties. 

The systematic uncertainties on the cross section ratio arise from five main sources: the neutrino flux (which 
largely cancels in the ratio), the neutrino interaction cross sections (which affect the background predictions), the 
target nucleon momentum distribution (which accounts for the model dependence of our unfolded neutrino energy), 
hadron re-interactions in the detector, and the detector simulation (which describes light propagation in the oil). In the 
region of highest statistics (about 1 GeV), there is roughly an 8% fractional error on the ratio resulting from hadron 
re-scattering in the detector, 6% from neutrino cross sections, 4% from the detector simulation, 2% from the nucleon 
momentum distribution, 2% from the neutrino flux, and 2% from the statistics of the two samples. 

The cross section ratio reported by all prior experimental measurements [7, 8, 9] has been one in which the effects 
of final state interactions (FSI) in the target nucleus have been removed using MC or deuterium was used as a target. 
Solely for the purpose of comparison, we now extract a similarly corrected value. The FSI-corrected ratio is defined 
as the ratio of CC\k^ to CCQE events at the initial vertex and before any hadronic re-interactions. Thus, the signal 
fractions and cut efficiencies for the FSI-corrected ratio include corrections for intra-nuclear hadron re-scattering based 
on the MC's model for nuclear effects [12]. The measurement proceeds exactly as for the observed ratio, except that 
now we define CCl tt"*" and CCQE, rather than CCl TT^-Iike and CCQE-like, events as signal for the respective samples. 

Here we limit our comparison to those experiments which reported both CCQE and CCl 7r+ cross sections, using the 
same energy bins for each of these interactions, so as to facihtate comparison with our measured CCln^ /CCQE ratio. 
Our result agrees with both ANL, which used a deuterium target, and K2K, which used C^H%. In order to make this 
comparison, the MiniBooNE and K2K results have been re-scaled to an isoscalar target. To perform this correction, we 
rescale the ratio by a factor of (1 — r)sp, where r is the ratio of neutrons to protons in the target and Sp is the fraction of 
7r+ production that is predicted (by MC) to occur on protons. The resulting scaling factor is 0.80 for MiniBooNE; for 
K2K we use the factor of 0.89 provided in [9]. The results have not been corrected for their differing nuclear targets 
nor for the application of explicit invariant mass requirements (although the latter are similar). ANL used an explicit 
cut on invariant mass W < 1.4 GeV [7] . While no invariant mass cut is used in this analysis, the MiniBooNE spectrum 
is such that CC \n'^ events occur only in the region 1^ < 1 .6 GeV; similarly, K2K's measurement covers the region W 
< 2 GeV [9J . 

The dominant reason for the difference between the ratios presented in Figure 1 is intra-nuclear pion absorption in 
CCl7r+ events, which cause these events to look CCQE-like. As a result of 7r+ absorption, a significant number of 
CCIk^ events migrate from CCtt"*" sample to CCQE-Vke, one. Thus, the FSI-corrected ratio, is 15% to 30% higher 
than the observed ratio in our energy range. 

In summary, MiniBooNE has measured the ratio of CCl7r+-hke to CCQE-\ik& events for neutrinos with energy 
0.4Gey <Ey < 2AGeV incident on CH2. This is the first time such a ratio has been reported. Additionally, the ratio 
of the CCl 7r+ and CCQE cross sections at the vertex has been extracted using MC to remove the effects of final state 
interactions, in order to facilitate comparison with previous experimental measurements. The results are summarized 
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FIGURE 1. Left: Observed CClTf^ -like/CCQE-\ike cross section ratio on CH2, including both statistical and systematic uncer- 
tainties, compared with the MC prediction [12], The data have not been corrected for hadronic re-interactions. Right: FSI-corrected 
CCln^ to CCQE cross section ratio on CH2 compared with results from ANL (D2) [7] and K2K (Cg/Zg) [9]. The data have been 
corrected for final state interactions and re-scaled for an isoscalar target . 



in Table I. The measured ratios agree with prediction [12, 13, 15] and previous data [7, 8, 9]. 
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